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Abstract. We study, using direct orbit integrations, the kinematic response of the outer stellar disk to the presence 
of a central bar, as in the Milky- Way. We find that the bar's outer Lindblad resonance (OLR) causes significant 
perturbations of the velocity moments. With increasing velocity dispersion, the radius of these perturbations is 
shifted outwards, beyond the nominal position of the OLR, but also the disk becomes less responsive. If we follow 
Dehnen (2000) in assuming that the OLR occurs just inside the Solar circle and that the Sun lags the bar major 
axis by ~ 20°, we find (1) no significant radial motion of the local standard of rest (LSR), (2) a vertex deviation 
of ~ 10° and (3) a lower ratio cri/ci of the principal components of the velocity-dispersion tensor than for an 
unperturbed disk. All of these are actually consistent with the observations of the Solar-neighbourhood kinematics. 
Thus it seems that at least the lowest-order deviations of the local-disk kinematics from simple expectations based 
on axisymmetric equilibrium can be attributed entirely to the influence of the Galactic bar. 
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1. Introduction 

Several observations since the 1980s have shown beyond 
any doubt that the Milky Way is actually a barred galaxy. 
Although there is strong evidence that bars in galaxies 
are confined to regions within the radius i?cR of the co- 
rotation resonance (CR), a bar may nonetheless influence 
the outer parts of its host galaxy, most obviously by res- 
onant phenomena. It has been shown by Dehnen (2000) 
and Fux (2001) that if the Sun is just outside the OLR 
such influence can explain the so-called u-anomaly. The 
latter is a bi-modality in the local velocity distribution of 
old stars, that has been inferred from HIPPARCOS data 
(Dehnen, 1998). 

The galactic bar might also have some relevance in ex- 
plaining the vertex deviation. This is the mis-alignment 
of the local velocity dispersion ellipsoid with respect to 
the local radial direction in the Galaxy and vanishes for 
any axisymmetric stellar dynamical equilibrium. Using 
HIPPARCOS data, the vertex deviation has been shown 
(Dehnen & Binney, 1998) to reach as high as 30° for young 
stellar populations, and to be ~ 10° for the old ones. 
Whereas the high vertex deviation of the young popu- 
lations is most probably due to moving groups, i.e. devi- 
ations from dynamical equilibrium, and thus in a sense 
accidental (Binney & Merrifield, 1998), there should be 
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dynamical reasons for its occurrence with the old popula- 
tions, i.e. deviations from axisymmetry. 

Another anomaly in the stellar kinematics observed 
in the Solar neighbourhood that may also be related to 
the bar's influence is the low value for the ratio ctI/ctJ of 
the principal components of the velocity-dispersion tensor 
of only 0.42j]Q;Q4 for old stellar populations (Dehnen & 
Binney, 1998). For a flat rotation curve for the Milky Way 
and in the limit of vanishing velocity dispersion, epicycle 
theory gives a lower limit of 0.5 for this ratio and the 
actual expectations from axisymmetric models are « 0.6 
(Evans & CoUett, 1993; Dehnen, 1999), significantly in- 
consistent with the data. 

In the present work, we will follow the approach of 
Dehnen (2000) further by investigating the velocity distri- 
bution in its spatial variability. Our primary interest is in 
the low-order moments of this distribution, i.e. the mean 
(streaming) velocity and velocity dispersion. In particular, 
we want to quantify whether the influence of the Galactic 
bar may explain the aforementioned anomalies, the vertex 
deviation and velocity-dispersion axis ratio, observed for 
the old stellar populations. Our approach applies to the 
kinematics in the solar neighbourhood, while at the same 
time it constitutes a completely general analysis of the bar 
influence in a stellar disk. 

The outline of the paper is as follows: In Sect. 2 we will 
describe our simulations, Sect. 3 gives a rather technical 
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presentation of the results, and Sect. 4 concludes the paper 
by a discussion of the results. 

2. Simulation of Bar Influence 

Since we do not want to construct a self-consistent model 
of the Galaxy, we just study the stellar dynamics for a 
simple model potential. In order to arrive at (stationary) 
equilibrium, we slowly add to the underlying axisymmetric 
Galactic potential the non-axisymmetric component of the 
bar (the bar monopolc is assumed to be already accounted 
for by the Galactic potential). We do not pay very much 
attention to the inner parts (inside co-rotation), and we 
neglect influences of vortical motion, so our model is two- 
dimensional. The calculation consists of orbit integration 
of a large ensemble of phase-space points representing the 
initial equilibrium. This numerical technique, which may 
be called restricted N-body method, is equivalent to first- 
order perturbation theory, since the self-gravity due to 
the wake induced by the perturbation (bar) is neglected. 
After the orbit integrations, the velocity moments at cer- 
tain times are computed from the phase-space positions 
of the trajectories. 




Fig. 1. Geometry of the rotation of the galaxy and definition 
of coordinate system. Note that positive values of radial ve- 
locity u are taken to point inwards and that azimuth angle tp 
is measured from the bar axis in the mathematically positive 
sense, but against the direction of bar rotatiou (modulo 180°). 
Also shown is a velocity dispersion ellipsoid with its principal 
components cri and cr2 and a (positive) vertex deviation Iv . 



2.1. Sampling 

The sampling of the initial phase-space points is done as 
described in Dehnen (1999), using the distribution func- 
tion (eq. (10) in Dehnen, 1999). 



f{E,L) 



n{RE)^{RE) 
ttk{Re)(t^{Re) 



exp 



Vt{RE)[L - L,{Re) 



'■{Re) 



(1) 

where fl{R), K,{R), and Lc{R) are the azimuthal and epicy- 
cle frequency and angular momentum of the circular orbit 
at radius R, while Re is the radius of the circular orbit 
with energy E. With this choice for the distribution func- 
tion, the coUisionless Boltzmann equation is satisfied at 
t = (since / depends only on the integrals of motion E 
and L) and the surface density and radial velocity disper- 
sion of the disk follow approximately (for a quantitative 
comparison, see Dehnen, 1999, and also Fig. 3) those given 
with the parameter functions E(i?) and cr{R), respectively. 
Here, we assiimc exponentials for both of them: 

S(i?) = Eoe^-^^-f^)/-"^, (j(i?) = troe^-f*"-")/-^". (2) 

If not stated otherwise, we choose R^ = Ro and i?s = 
0.33i?o by standard, Rq being the distance of the Sun from 
the galactic center. If not stated otherwise, ao = 0.2uo, 
where vq is the circular velocity at Rq. In this way, samples 
of K = 10^ initial phase-space points are created. 



2.2. Model Potential and Orbit Integration 

Orbit integration and adiabatic growth of a quadrupole 
bar is done similarly to Dehnen (2000). The galactic back- 
ground potential is chosen to give a power law in the ve- 
locity curve 

v,iR) = vo [R/Rof , (3) 



for /? / 



(4) 



namely: 

[ In (R/Ro) for /3 = 

For the bar potential i>i = ^i{R,(p,t), we only use a 
quadrupole, since higher poles are much less important 
at large radii, following Dehnen (2000): 



f 2 - (R/Rb) for R<Rb 
$1 = -Ab{t) cos2ip X < 

[ (Rb/R) for R> Rb 



(5) 



where Rb is the size of the bar. The angle (p is defined in 
the frame rotating at pattern speed Qb (see Fig. 1). The 
strength Ab of the bar is increased from to a value Af 
during time < t < ti according to 



5^3 1 
-£ H S. + - 



,1 



1 



(6) 

and stays constant a.t Ab = Af after ti. With this func- 
tional form, $1 and its first and second time derivative are 
continuous, thus representing adiabatic growth of the bar 
and ensuring a smooth transition. The final bar strength 
Af is controlled via the dimensionless model parameter 
ao = C({Ro) with 

A / n \ 3 

a{R) = 3^ 



Rb 



(7) 



which is the ratio of the forces due to $i and i>o at galacto- 
centric radius R > Rb on the bar's major axis. Our stan- 
dard choice is a.^ = 0.01, which is consistent with the data 
for the Milky Way. 

Orbit integration is performed in the co-rotating frame 
up to a time t2 > ti using a 5th order integrator. In con- 
trast to Dehnen (2000), we integrate forward in time. 
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2.3. Calculation of Velocity Moments 

The {p, q)th moment of the velocity distribution at posi- 
tion r is defined as 



Mpq{r,t) ^ J dudv u^v"^ f{t, x = r, v) (8a) 
= J dxdv S{r - x) u^v'^ f{t, X, v) (8b) 

with V = {u,v), where u and v denote, respectively,the 
radial and azimuthal velocity component (see Fig. 1). We 
do not know the value f{t, x, v) of the distribution func- 
tion at any time t > 0, but instead have a representa- 
tive sample {xk{t),Vk{t)} of phase-space points. Hence, 
we compute the moment integral (8b) via Monte-Carlo 
integration, resulting in 



Mpg{r,t) 



tot \ ^ 



1 



K 



^4 



\Xk 



fife 



(9) 



where we have replaced the 5 function in Eq. (8b) by 



\xk — ■'"l/cfc) with the weight function 



w{d) = -(1 -d^) 0(1 -d-^), 

TT 



(10) 



where 9 is the Heaviside function. The parameter gives 
the radius over which the fcth trajectory contributes to the 
moment integrals. We adjusted tk such that a constant 
number of sampled orbits was expected to fall in the area 
of radius centered on Xk by the assumed exponential 
surface brightness distribution of the disk. 

In practice, the zeroth, first and second moments are 
estimated in this way and the mean velocities and ve- 
locity dispersion tensor^ are then evaluated from fipq = 
Mpq/MoQ as M = /iio, w = Hoi and cr^„ = ^20 - A^?o> 
^Iv = /^02 - Moi' '^uv = Mil - MioMoi- 

Once the orbit integration is done, one can easily 
switch to a model with initial distribution function fi^f 
by weighting each orbit with the ratio fi/f, which ac- 
counts for the fact that the trajectories were actually 
sampled from /. However, for this method to be useful, 
the ratio /i// must not become too large, because oth- 
erwise the moment estimates are dominated by a few or- 
bits with large weights. Since, by virtue of the collision- 
less Boltzmann equation, the values of the distribution 
functions are conserved along trajectories, the ratio fi/ f 
is conserved, too, and can be evaluated at time t = 0, 
for which the distribution functions can be computed via 
Eq. (1). In this way, we did a switch to exponential disks 
with velocity dispersions ctq < 0.2vq (the value of the sam- 
pling distribution function), such that /i// < [0.2'i;o/o'o]^ 
(this follows from Eq. (1); for ctq > 0.2?;o, the ratio /i// 
can become very large) 



^ Notation: a'^„ etc. are components of the tensor cr^, whose 
eigenvalues are the squares of the velocity dispersions in the 
principal directions. Note that a^^ may well be negative. 
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Fig. 2. Evolution towards stationarity as measured by the 
symmetry (13) for different models (arbitrary units, horizon- 
tal axis is time in bar rotation periods). Models shown have 
RolrJRo ~ 0.8 (squares) and Roi,r/Ro = 0.92 (triangles), 
and usual (filled) or long integration time (open). Bar growth 
is taking place in the first 5 periods for normal and in the 
first 10 periods for the long integration time. In addition, the 
value for the sampled starting distribution at t = is shown, 
indicating the noise level. 

2.4. Fourier Components of Velocity Moments 

In order to most easily follow the evolution of the disk, it is 
useful to do a Fourier transform on the azimuthal angle ip. 
We apply this to the mean velocities and dispersion tensor 
elements after we calculated these in the usual way (9), 
except that this time we do not use a weighting function, 
but take bins in R and (p to make sure every calculated 
trajectory contributes to the final result. We use a discrete 
Fourier transform of the following kind: 



N 



(11) 



i=i 



where ipj — jn /N . This gives an approximate Fourier ex- 
pansion 



JW ~ y + y cosnv?- 



n-l 



(cm COS mip + Sm sin mip), 
(12) 

where n = N/2 and N is supposed to be even. This has 
been done for iV = 32 in 25 radial bins. 

We also construct an estimate of the surface density E 
by dividing the number of sampled orbits per bin by the 
segment area of the bin, and apply the Fourier transform 
to this quantity as well. 
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2.5. Error Estimation 

We employ a bootstrap method: the calculations leading 
from the data set of every radial bin to its Fourier coeffi- 
cients arc redone for arbitrary subsamplcs of this data set. 
The rms scatter in the outcome of many such calculations 
gives a measure of the error. 

2.6. Symmetries and the Question of Stationarity 

In order to check, whether or how far our distributions 
have reached a stationary equilibrium state, we employ a 
symmetry consideration. As noted by Fux (2001), velocity 
distributions in an m = 2 symmetric potential have a 
symmetry 

f{R,ip,u,v) = f{R,TT - ip,-u,v), (13) 

which is a consequence of the azimuthal m = 2 symmetry 
and the time-reversal symmetry of stellar dynamics in con- 
junction with stationarity of the distribution function. We 
can turn this argument around and use the degree of sym- 
metry as a measure for stationarity. A quantitative mea- 
sure of the symmetry can be obtained in the obvious way 
by subtracting corresponding distribution function values 
over a grid in phase space and doing a quadratic sum over 
the grid points. Results of this kind of analysis are shown 
in Fig. 2. As can be seen, the simulations show only a 
small increase in asymmetry of 15% above the noise level 
(as given by the initial state as well as a very long-time 
simulation), and this is dropping rapidly after bar growth 
is finished. 

3. Results 

Kuijken & Tremaine (1991) gave an analytical expression 

for the behaviour of the mean velocities under the influ- 
ence of a non-axisymmetric perturbation of multipole or- 
der m in a linear approximation. For our case, this yields 

_ ^ ajR) SQb - O . 

3R (fi - niLR)(0 - Oolr) ^'"^ 

(14) 

a{R)v^ Aftb-n 

6R (S2 - S2iLR)(i2 - iioLR.) 

where Oilr and Oqlr refer to the circular frequencies at 
the radii of the inner and outer Lindblad resonance. For 
a flat rotation curve, = vq/R. In our range of interest, 
Eqs. (14) are dominated by a pole at the OLR. 

The results of our simulation, presented below, agree 
roughly with these expectations: the mean velocities show 
modulations u oc sm2(p, v oc cos 299. The sign of the mod- 
ulations varies with radius, sign changes should indicate 
resonances. For the dispersion tensor the situation is sim- 
ilar: the off-diagonal component shows a sin 2(fi perturba- 
tion, the diagonal components a cos 2ip. 



3.1. Fourier Analysis 

Since we have an azimuthal m = 2 symmetry, all coeflS- 
cients of odd index m are expected to vanish, what they 
do within their errors. Furthermore, in the undeveloped 
state before bar growth, all coefficients except those with 
m = are, by construction, zero within their errors. 

3.1.1. The m = Components: 

The dashed lines in Fig. 3 represent the initial axisym- 
metric case, which shows an exponential decline of S as 
well as in the diagonal components of velocity dispersion 
tensor. The mean azimuthal motion deviates from the cir- 
cular speed vq by the asymmetric drift, which is stronger 
for large ctq, as expected. The mean radial motion van- 
ishes, as required for any stationary model. 

The solid lines in Fig. 3 show the radial run of the 
m = components in the barred case. Apart from the 
components for u and c^^. which have to vanish for sta- 
tionary models, some signs of perturbation are visible, the 
stronger the smaller ctq. 

First of all, we learn that dispersive effects are quite 
efficient in drawing resonance-induced features away from 
the actual position of the resonance. Whereas for small 
velocity dispersion the association of the features with the 
OLR at i?oLR = 0.92i?o is clearly visible, they appear in 
the high dispersion case at a radial range of 1.1 to 1.4i?o, 
where naively one would not attribute them to the OLR. 

In particular, we have a bump in the v-cmve outside 
of OLR, which becomes more explicit with decreasing dis- 
persion while roughly keeping its absolute magnitude. In 
contrast to this, the asymmetric drift as the dominant de- 
viation of V from the nominal rotation velocity vq dimin- 
ishes with decreasing dispersion. There is a similar bump 
in (T^^, whereas we have a two- fold feature in a^y. 

Most of the perturbative features in the m = com- 
ponents are induced by the opposite orientation of the 
near-circular orbits on either side of the resonance. 

3.1.2. The m = 2 Components 

Because of the symmetry of the problem, we expect the 
m = 2 components to be the dominant ones. As for every 
m > 0, we have one more degree of freedom here, since 
we have a cosine- and a sine-component (see Fig. 4), or 
equivalently, an amplitude and a phase. 

The cosine-terms of u and cr^^ as well as the sine terms 
of S, V, (7^„, and a^^ vanish (within their error), which is 
to be expected from Eq. (13) for a stationary model. The 
behaviour of the remaining non-vanishing terms agrees 
roughly with the expectation from simple linear theory, 
cf. Eq. (14). However, there are two significant deviations. 
First, as ctq increases (from left to right panels in Fig. 4), 
the resonant features are shifted away from i?oLR to larger 
radii and are also somewhat smoothed out. Second, there 
is an additional feature at ii ~ O.6R0, in particular, for 
V. When inspecting the orbits of stars dominating at this 
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Fig. 3. Axisymmetric (m = 0) components of velocity distribution in the undeveloped state before bar growth (dashed) and 
in the final state (solid), plotted against radius. Results for models with initial ao — O.OSno, O.lvo and 0.2vo are shown in the 
left, middle and right panels, respectively. Velocities are given in units of vq and E in arbitrary units. The model shown has an 
OLR- radius of 0.92iio, corresponding to a CR of 0.539i?o (dashed vertical lines). 



radius, we find that many belong to an orbit family asso- 
ciated with the stable Lagrange points L4 and L5, see Fig. 
5. These orbits can reach radii far beyond co-rotation and, 
in a frame co-rotating with the bar, perform a retrograde 
motion around the Lagrange points. They thus result in a 
reduced v at azimuths perpendicular to the bar and hence 
lead to a positive cos2(^ component. 

3.1.3. Higher-Order Components 

The m — 4 components are usually smaller by at least 
a factor 3 compared to the m = 2 modes, but they are 



significantly different from zero. In contrast, the analytical 
model (14) had no excitation of higher modes at all, since 
mode coupling is not included in a linear approximation. 

The excitation of the m = A modes follows the overall 
pattern seen in the m = 2 case: in accordance with the 
symmetry requirements, we have a sin(4<y9) in u and cr^^, 
and a cos(2(/3) in all the other components. Our signal-to- 
noise ratio is in the range of about 3 for the m — A case, 
and less for higher components, so surely the significance 
of our model is petering out. Nevertheless, some m = 8 
modes do seem to be discernible still. 
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Fig. 4. m = 2 Fourier cosine (solid) and sine (dashed) components, plotted against radius. Panels as in Fig. 3. For surface 
density E, the value shown is relative to the undisturbed m — value. OLR and CR are at 0.92i?o and 0.539i?o respectively 
(dashed vertical lines). 



3.2. Variation of Parameters 

Doubling the bar strength leads to a doubling of the mag- 
nitude of the effects. This is true for the differences be- 
tween the undeveloped and final state as well as for Fourier 
coefficients. 

Changing the bar size, i.e. the parameter i?t, in (5), 
produces some local variations in the magnitude of the 
effects, but it does not shift them in location. Of course, 
this is to be expected if the location is determined by 
resonance conditions. 

Variation of the disk scale length does not seem to 
have any strong effects, which is to be expected, since it 



only slightly affects the sampling of the orbits but not the 
dynamics. 



The shape of the velocity curve has a mild effect: we 
tried models with a slightly rising (/3 = O.f) or falling 
((3 — —O.f) velocity profile, and found very similar results 
as for /3 — 0. (Note that we kept the OLR radius fixed, so 
models with different velocity profiles have different rota- 
tion frequencies of the bar.) This is expected, too, since 
(3 essentially controls the distances between various res- 
onances (growing wider for larger /3), while most of our 
results are dominated by a single resonance, the OLR. 
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Fig. 5. Orbit near co- rotation that circulates the L4 point in a 
retrograde sense (anti-clockwise for a clockwise rotating bar). 
The bar is indicated by a line, while the dotted circles corre- 
spond to CR and OLR. In our model, a considerable fraction 
of stars is trapped in orbits of this family, resulting in the 
somewhat peculiar run of the velocity moments (Fig. 4). 

3.3. The Axis Ratio of the Velocity Ellipsoid 

Fig. 6 shows the ratio of the principal axes of the velocity 
dispersion ellipsoid, i.e. the eigenvalues 2 of the tensor 
cr^. For a flat rotation curve, this ratio is often expected 
to be 0.5. Indeed, in the undisturbed case (for which cti = 
auu and 0-2 — cjyy) 

lim ^ = — ^ = - 1 -I- , (15 

where k is the epicycle frequency. This relation is often 
given without the limit and then referred to as Oort's 
relation. It is, however, important to note that the error 
in Oort's relation is considerable already for cto — 0.2wo 
(Evans & CoUett, 1993; Dehnen, 1999), which means that 
using it for the old stellar disk of the Milky Way is, at best, 
dangerous. This can be seen clearly from Fig. 6 showing 
that cr| / a\ for the unperturbed case (dotted) significantly 
deviates from 0.5 for a warm stellar disk (cto = 0.2i;o). 

After formation of the bar perturbation, we have large 
variations in this quantity. For ctq = 0.2tio and at i? = 
i?o, the values are generally somewhat higher than Oort's 
value 0.5 but smaller than for an undisturbed disk. For 
directions that are roughly along the bar [^p ~ 0), this 
ratio rises sharply just outside the solar circle to values 
reaching as high as 0.8. 

It is instructive to take a look at the same quantities 
for a smaller velocity dispersion of ctq ~ 0.05wo or O.Iwq 
(left and middle panels of Fig. 6). First of all, the axis ratio 
in the undisturbed case (dashed line in the right panels) 
is much closer to Oort's value of 0.5 here, as expected. 
Generally, the observed features have smaller width, i.e. 



they are less washed out by dispersion, but the effect of the 
bar, i.e. the difference between the solid and dashed lines, 
is still even quantitatively similar for the different veloc- 
ity dispersions. Thus, the deviation of o\la\ from Oort's 
value may be decomposed into a velocity-dispersion de- 
pendent term, which consists of a general elevation of only 
mild radial dependence and a bar-induced term, which 
varies spatially and appears to be negative for the Solar 
position. 

As a consequence, for mildly warm stellar disks (cr ^ 

0. 2uo), the ratio (j\lu\ in the Solar neighbourhood may 
well drop below Oort's value. 

3.4. Vertex Deviation 

The so-called vertex deviation 

1 2cr2 
^„ = - arctan ^ ^ (16) 

Z fT — fT 

uu vv 

is the angle between the direction of the largest velocity 
dispersion and the line to the Galactic center, see also Fig. 

1. For axisymmetric equilibrium models, £.v — 0. 

In Fig. 7, we plot £„, computed from Fourier coeffi- 
cients up to m = 2 only in order to eliminate short-scale 
fluctuations. We find that the bar-induced vertex devia- 
tion decreases with increasing velocity dispersion from up 
to - 30° for (To = O.Obvo to 5° - 10° for ctq = 0.2vo. In 
the direction of the bar and perpendicular to it, £v = 0, 
as expected from symmetry. For moderately warm stellar 
disks ((To ~ O.lwo) the vertex deviation at azimuth angles 
between 0° and 90° (including the Solar azimuth) is pos- 
itive at and inside the solar circle and negative for some 
range outside. The actual radius where the sign change 
occurs is obviously coupled to the OLR, but again shifted 
to the outside. For a cold stellar disk, the bar-induced ver- 
tex deviation displays are more complicated pattern and 
may reach as high as 40°. 

Obviously, the vertex deviation is antisymmetric with 
respect to if 1— s- (tt — tp), whereas the velocity disper- 
sion axis ratio is symmetric. This is a simple conse- 
quence of symmetry (13), which in particular implies 

/2 2 2^ /'2 2 2\ 

\'-^ uu T '-^ uv ' ^ vv ) ' * l^uu ' ~^m> ' "'iidJ' 

4. Discussion 

We have studied the bar-induced variations of the mean 
velocities and velocity dispersions in the outer disk of a 
barred galaxy with special emphasis on the situation in 
the Milky Way. We generally find that these variations 
are largest near the outer Lindblad resonance (OLR), the 
apparent radius of which, however, may be shifted out- 
wards by 10-20% due to the non-circularity of the orbits 
in a warm stellar disk. 

Radial motions u of standards of rest can be seen to 
occur quite frequently (see Fig. 8), and can reach magni- 
tudes of the order of about 0.02wo, corresponding to about 
5 kms^^ for the Milky Way. Because of the sin2(/3 depen- 
dence, these would be maximal at ip = 45°, which is quite 
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ratio al/o"! of the principal components of the velocity-dispersion tensor for the final [solid) and the initial 
state (dashed) at various azimuths for the models with Rqlb. = 0.927?o and ao of 0.05i;o (lejf), O.lvo (middle), and 



near the proposed position of the Sun of « 30°. In its 
radial dependence however, u swings through zero shortly 
outside of the OLR, and it may well be that the Sun just 
meets that point. So we cannot give a definite prediction 
for the bar-induced radial motion of the LSR here, not 
even by sign only, except that it should be very small (at 
most a few kms~^). 

Variations of the mean azimuthal velocity are also 
present, but here we have to deal with several perturbation 
effects, the dominant one for realistic cases being tradi- 
tional axisymmetric asymmetric drift. We can construct a 
picture of the bar-induced t)-modifications only (see Fig. 



8) by subtracting v in the undeveloped case before bar 
growth, thus also canceling asymmetric drift. However, 
this is, of course, impossible for real galaxies, where only 
the combination of circular speed (LSR motion), asym- 
metric drift, and bar induced drift is measurable (in prin- 
ciple) . In observational data, the effect of asymmetric drift 
can be corrected for, as done by Dehnen & Binney (1998), 
by calculating v for stellar populations with different ve- 
locity dispersions and extrapolating to vanishing disper- 
sion. This approach, however, cannot be applied to the 
bar- induced perturbations in v, since, as can be seen in 
Fig. 3, the relative magnitude of the bar- induced wiggle 
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deviation iy at various azimuths for the same models as in Fig. 6. The highly uncertain outliers at ~ 1.2Ro, 
= O.lvo are caused by a near-circular velocity ellipsoid (cri ~ cr2, see Fig. 6), which renders the vertex deviation 



keeps constant at about 2%, independent of velocity dis- 
persion. So, as long as only local measurements are avail- 
able, there is no way to discern the bar-induced azimuthal 
velocity perturbations from variations in the background 
rotation curve, and so it will be rather hard to draw any 
observational constraints from the mean azimuthal veloc- 
ities. 

For an axisymmetric equilibrium, the no vertex devi- 
ation, i.e. iy = 0, is expected, while in the Solar neigh- 
bourhood £v is observed to drop from £y w 30° for blue 
to w 10° for red main-sequence stars (Dehnen & Binney, 
1998). The traditional interpretation of these vertex devi- 
ations is that young stars deviate from equilibrium, and 



thus even in an axisymmetric galaxy give rise to £y ^ 0, 
which due to the decreasing number of young stars at later 
stellar types is a decreasing function of stellar colour. 

Our calculations show that for an equilibrium in a 
barred Milky Way, with bar orientation, strength and pat- 
tern speed consistent with other data, a vertex deviation 
of the size and direction as observed emerges naturally. 
Moreover, we also found that for dynamically cooler sub- 
populations, i.e. bluer or younger stars, the bar-induced 
vertex deviation increases in amplitude, very similar to 
the observed values. 

This gives strong support for the hypothesis that the 
vertex deviation observed in the Solar neighbourhood is 
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Fig. 8. Bar-induced deviations of the mean velocity from the unperturbed state up to m = 2 for ao = 0.2vo. The bar is aligned 
with the horizontal axis and is supposed to rotate clockwise. The solar circle is dashed, dotted circles are CR and OLR. 



predominantly caused by deviations from axisymmetry 
rather than from equilibrium. This explanation also natu- 
rally accounts for the fact that £y for young stars has the 
same direction as for old ones, which with the traditional 
explanation would be a chance coincidence. 

The axis ratio of the (principal components of the) 
velocity dispersion tensor, crl/u^, is clearly affected by the 
central bar. In particular, values less than 0.5, Oort's value 
for a flat rotation curve, are possible (Oort's value is a 
lower limit for an axisymmetric galaxy, see Sect. 3.3 and 
Evans & CoUett, 1993; Dehnen, 1999). This nicely fits to 
the values inferred from HIPPARCOS data (Dehnen & 
Binney, 1998), which give frl/crj ^ 0.42 for the old stellar 
disk. 

In the Milky Way, as well as in many barred galaxies, 
the bar is not the only deviation from a smooth axisym- 
metric background. Most prominently, spiral arm struc- 
ture and an elliptic (or oval) disk and/or halo add further 
non-axisymmetric perturbations. However, we will defer 
the discussion of the influence of spiral arms on the outer 
disk kinematics to a future paper. 
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